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We investigate a dilute Fermi gas of polar molecules confined into a bilayer setup with dipole 
moments polarized perpendicular to the layers. In particular, we consider the extreme case of 
population imbalance, where we have only one particle in one layer and many particles in the 
other one. The single molecule is attracted by the dilute Fermi-gas through the inter-layer dipole- 
dipole force presenting an interesting impurity problem with longrange anistropic interaction. We 
calulate the chemical potential of the impurity, in second order perturbation theory and in ladder 
approximation with a Brueckner-Hartree-Fock approach. Moreover, we determine the momentum 
relaxation rate of the impurity, which is related to the "dipolar" drag effect. For a confined system 
we relate the results for the chemical potential with the measurement of the collective modes of the 
impurity. The momentum relaxation rate provide instead an estimate on how quickly the oscillations 
are damped. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.-d 



I. INTRODUCTION 



During the last years ultracold dipolar gases have at- 
tracted great interest because the dipole-dipole interac- 
tion (DDI) drastically changes the nature of quantum de- 
generate regimes compared to ordinary short-range inter- 
acting gases [H, ■ In recent experiments polar molecules 
in the ground ro- vibrational state has been created [!, H| 
and cooled towards quantum degeneracy Q. The main 
obstacle is the decay of the system due to ultracold chem- 
ical reactions found in recent experiments However, 
if the molecules are confined in a 2D geometry and ori- 
ented perpendicularly to the plane of their motion by a 
strong external electric field, these reactions are expected 
to be suppressed by the intermolecular repulsion [6[ . Bi- 
layer arrangements with this dipolar orientation are par- 
ticularly interesting, since they allow for both, stabil- 
ity against chemical reactions and effects induced by the 
anisotropy of the DDI. The attractive inter-layer inter- 
action gives rise to a peculiar two dimensional scattering 
behaviour and bound states JtHSJ and may induce inter- 
esting phenomena like inter-layer pairing and superflu- 
idity 04H| and non-local state-changing collisions [l^ . 
Notice that the bilayer system can be also thought as 
a single layer system with two different species interact- 
ing through the intra-layer and the inter-layer DDI. Until 
now only the balanced situation with an equal number 
of particles in both layers has been studied. 

For atomic short-range interacting Fermi gases of dif- 
ferent species (e.g. atoms in two different hyperfine 
states) however, imbalanced and in particular highly im- 
balanced gases have been extensively studied experimen- 
tally [l3[ and theoretically (see, e.g., the recent review 
[lH and references therein). The building block for the 
understanding of such systems is the solution of the lim- 
iting case of a single impurity atom interacting via a 
short range potential with an ideal atomic Fermi gas. 
Such a problem is not only relevant in the field of ultra- 
cold gases, since it is related to the so-called impurity 




FIG. 1: Scheme of the system under consideration. 



problem, which is present also in other areas of physics. 
Analogously it is interesting to study the bilayer system 
of polar molecules with a high population imbalance, re- 
sembling an impurity interacting via the inter-layer DDI 
with a Fermi gas (see Fig. [1]). In the following we call 
the single particle in one layer impurity or Fermi-polaron 
or polaron in analogy with electrons in a crystal dressed 
by the bosonic (phonon) bath. 

An important quantity characterizing the impurity is 
its chemical potential, namely the (in our case negative) 
energy difference of the ground state with and without 
the impurity at rest. We will sometimes call this quan- 
tity also interaction energy or binding energy. For short- 
range interacting Fermi gases this quantity has been suc- 
cessfully calculated with various methods in one [HI, [l6[ , 
two (i~7H20| and three dimensions [2l| - [23| . 

In this paper we investigate the problem of an impu- 
rity in one layer (A) interacting via the DDI with the 
Fermi gas phase in the other layer (B) as sketched in 
Fig. [I}. Assuming no inter-layer tunnelling we calculate 
the chemical potential of the impurity in second order 
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perturbation theory and in ladder approximation using 
a Brueckner Hartree Fock approach. It has been shown 
in [21], H3 that for short range interacting atomic gases 
the latter is equivalent to a self-consistent T-matrix ap- 
proach and a variational ansatz, which have been applied 
succesfully to the polaron problem. Moreover, we study 
some dynamical properties of the impurity: we estimate 
the momentum relaxation rate and calculate the eigen- 
frequency of the impurity in an external trap. Note that 
we consider through the whole paper weak interactions 
and a dilute Fermi gas of polar molecules. 

The structure of the paper is as follows: in section Hi Al 
we introduce the relevant physical quantities and study 
the polaron energy in second order perturbation theory. 
In section IIIBI we extend our analysis to higher orders 
and calculate the chemical potential in ladder approxi- 
mation within Brueckner-Hartree-Fock approach. In sec- 
tion III we estimate the momentum relaxation time of a 
moving impurity. Finally in section IV we use our re- 
sults for homogeneous systems to investigate the motion 
of the impurity in an external trap using the local density 
approximation. 



II. CHEMICAL POTENTIAL 



A. Second order perturbation theory 



whose Fourier transform is 



V(q) = V AU {<\) = -2nd 2 qe 



-Xq 



(3) 



In first order perturbation theory the interaction en- 
ergy reads = (GS\H AB \GS), where \GS) = 
{"pi.s} , P) is the ground state of the system with Nb 
particles in layer B and one particle (the impurity) in 
layer A with momentum P. The first order contribution 



£?« = n B V(q = 0) = 



(4) 



where ns = Nb /V is the particle density in layer B, van- 
ishes. Indeed the potential is partitially attractive and 
partially repulsive such that it satisfies f drV (r) = 0. 
Similiarly, the second order contribution can be found as 



pi 



q#0 11 



2|^(q)| 2 (l-n pi+q , B ) 
f P 2 -( P i+q) 2 -(P-q) 2 ' 

(5) 



where rtk,s = ©(^f.b - k) denotes the Fermi-distribution 
at zero temperature. At rest (P = 0) this energy is the 
chemical potential of the impurity and in Sect. Ill B[ we 
discuss it as a function of the dipole moment d, the inter- 
layer distance A and the Fermi-momentum kF,B of the gas 
in layer B. 



The model Hamiltonian can be written as 



B. Brueckner-Hartree-Fock approach 



pi,p 2 ,q,<T,a-' 



(1) 



pi +q,cr C p 2 -q,cr' C P2 C Pl >o- ' 



where tk is the single particle energy and c pcr (c pcr ) cre- 
ates (annihilates) a fermionic particle with momentum 
p in layer a = {A,B}. The inter-layer and intra-layer 
dipolar interaction potentials are V aa and V AB , respec- 
tively. Note that up to second order the inter-layer and 
the intra-layer interaction decouple, such that the en- 
ergy of the impurity does not depend on the intra-layer 
interaction. At higher orders, processes mediated by the 
intra-layer interaction might become important and have 
to be taken into account in the calculation of the chemical 
potential of the impurity. Thus in the following second 
order calculation we consider only the inter-layer inter- 
action and we can omit the indices of the potential. The 
intra-layer interaction energy has been calculated in first 
order perturbation theory in Refs. 0, [25[ 

For two molecules with mass m and dipole moment d 
placed in two layers at distance A the DDI potential reads 



V AB {r)=d z 



(2) 



(r 2 + A 2 ) 5/2 ' 

with r being the relative distance in the plane of motion, 



We calculate the polaron energy summing up the lad- 
der diagrams for the inter-layer interaction only, which 
at low enough density (in layer B) should give the most 
relevant contribution. The basic equation is the Bethe- 
Goldstone integral equation for the reaction matrix [27| . 
also called effective interaction (e.g., for the 2D electron 
gas [III). In our case the Bethe-Goldstone integral equa- 
tion for the effective interaction between a particle in 
layer B with momentum ki and the impurity atom in 
layer A with momentum k2 can be written as 

f dk 

<7(e p ,k 1) k 2 ;q) = ^(q) + J _-_y(| q - k|) x 

;(e p ,ki,k 2 ;k) 



2m (1 - n kl+k , B ) 

U2 t L2 



2me p + kl + kl- (ki + k) 2 - (k 2 - k) 



(6) 



In Eq. © q is the momentum transfer, V(q) is the 
Fourier transform of the inter-layer potential ((31), the 
Fermi distribution function at zero temperature. Note, 
that the interaction energy e p has been included in the 

k 2 k 2 

initial energy of excitation processes e p + w 1 - + • This 
interaction energy or correlation energy follows then self- 
consistcntly from the mean value of the effective interac- 
tion e p = (g(e p ,ki,k 2 ;q)) 

For k 2 = one gets the rest correlation energy e° of 
the polaron, and by expanding this solution in k\ one 
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get its effective mass m* as usual by the relation E — 
e p + k%/2m*. We remind that in Eq. © for the effective 
interaction only ladder diagrams are summed and the 
Fermi sea limits the momenta in the intermediate states. 

For shortrange interacting Fermi gases the integral 
equation ^ cannot be solved directly due to the infi- 
nite hard core of the potential. Especially in 2D one has 
to renormalize the equation by replacing the potential 
through the off-shell scattering amplitude and solve the 
renormalized equation. We followed this strategy in Ref. 
)_19f. However, the function of the inter-layer DDI V(q) 
as given by Eq.Q goes to zero for small as well as for 
large arguments, such that it does not cause any diver- 
gencies. This allows for a direct numerical solution of 
Eq. ([6]) in our case. Then we can compare this solution 
with the one obtained by means of the renormalization 
strategy as discussed in |19| . 

For our problem the Bethe-Goldstone equation © is a 
two-dimensional integral equation in the last argument of 
g(e p , ki, k2; q) which for &2 = can be written in units of 
kp t B as a Fredholm integral equation of the second kind 

g(elh;q,6) = V(q)+ 



dk r 2 * dcj> 
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with the angles 9 = Z(ki,q) and 4> = Z(ki,k). In 
Eq. Q we have introduced the dimensionless inter-layer 
dipolar potential 



V (q) = -2ir(r*k F , B )e 



-(\k FiB )q 



(8) 



with r* = d 2 m/H 2 the length scale associated with the 
dipolar interaction and the kernel 



K(e p ,ki;q,9,k,< 



fc^(|q-k|){l-0(l-|ki+k|)} 



e ,/ {2cf,b) — k 2 — kk\ cos < 



(9) 



The mean value of the solution of Eq.Q gives the 
chemical potential e° = e^/{2eF,B) of the impurity self- 
consistently 



e° = 
p 



1 dkiki f 00 dqq 



27r d6 

(2tt) Jo (2ir)J ~(2w) 
l-e(l-|k 1+ q|) 

£p ~ I 2 ~~ 9^1 cos 4> 



V(q). 



(10) 



g(i p ,k i;q ,e). 



We can take as quantities characterizing the inter-layer 
DDI the dipolar strength r*kp,B and the inter-layer dis- 
tance \kp t B- The quantity r*kF,B is also the strength of 
the intra-layer interaction in layer B for which we con- 
sider only small values, such that we are far away from 
the liquid- crystal phase transition [2!|. In this case the 
intra-layer interaction could be taken into account by in- 
troducing renormalized Fermi liquid Landau parameters, 



which accounts, e.g., for an effective mass of the fermions 
111 layer B @. 

Fig. 2 shows the numerical solution in ladder approx- 
imation (Eqs. 0,© and (HUJ)) for r*k F , B = 0.5. The 
polaron energy is plotted as a function of the inverse 
inter-layer distance l/(Afci?.s)- For comparison the sec- 
ond order perturbation theorv result as obtained from 
Eq □ 
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FIG. 2: Main: polaron energy as a function of the inverse 
inter-layer distance for dipolar strength of r*kF,B = 0.5 cal- 
culated in second order perturbation theory (black dashed 
line) and in ladder approximation (red crosses). Inset: Zoom 
onto the weakly interacting regime, where second order per- 
turbation theory is a good approximation. 



As mentioned above one can also solve Eq. ^ approx- 
imately by replacing the potential through the scattering 
amplitude. In Ref. [19] we obtained in that way for the 
shortrange interacting Fermi gas an approximate analyt- 
ical solution for the polaron energy. The energy does not 
depend on the particular form of the interaction potential 
but only on the binding energy of the two-body bound 
state. We can apply the results of [l9| for all potentials 
providing a normal scattering behaviour, governed in 2D 
by a logarithmic on-shell scattering amplitude. Scatter- 
ing and bound states of the inter-layer dipolar potential 
has been analysed with different methods in detail in 
0, Q and one has that for a certain regime of parame- 
ters r* and A the low-energy scattering is indeed of the 
normal logarithmic form. In this case the polaron energy 
reads |19j 



-2e F , 



In 



(11) 



where \eb\ is now the binding energy of the two-body 
bound state of the inter-layer potential Eq. (J3J. This 
binding energy depends on r* and A and has been cal- 
culated in 0] and § as well. Using this results we can 
calculate the polaron energy from Eq. (TTT|) in the regime 
where the inter-layer potential provides a logarithmic low 
energy scattering amplitude. This condition is satisfied 
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FIG. 3: Polaron energy in log scale as a function of the inverse 
inter-layer distance for a dipolar strength of r*kp t B = 0.5 
calculated from Eq. [10] (red crosses) and from Eq. [TT] (black 
dashed line). The blue dashed dotted line indicates the two- 
body binding energy et for the inter-layer DDI potential (|3]) 
as calculated in Q 



for r*/X > 0.7 0. Fig. 3 illustrates the comparison 
of the results obtained from the numerical solution of 
Eqs. ©,© and (Unj) and formula JITJ for r*k F . B = 0.5 
and shows moreover the two-body binding energy for the 
inter-layer potential as obtained in Q. 

Since the self-consistent approach we use in this paper 
is equivalent to a variational ansatz with a polaron many 
body wave function, we expect that it cannot produce ac- 
curate results in the regime of strong two-particle interac- 
tion, where the particles are deeply bound in a molecular 
state. The investigation of this molecular regime is be- 
yond the scope of our paper. Moreover it has been shown 
for short-range interacting Fermi gases that the molec- 
ular regime in 2D, differently to the three-dimensional 
case, occurs for very large interaction strength [18l. |30|. 
For our case strong particle interaction means that for 
fixed r*kp^B the inverse inter-layer distance is very large. 
We find that for 1/(AA;f,s) ^> 1 the polaron energy con- 
verges to e° w | ej, | — 0.55eF.s, which is in good agreement 
with the result in [l8| for the polaron ansatz in the two- 
dimensional shortrange interacting gas. 



III. MOMENTUM RELAXATION RATE 

In the previous section we have calculated the energy 
of the polaron as a function of the interaction strength. 



We assumed the polaron is well defined with an energy 
given at finite momentum e q = e° + q 2 /(2m*), where 
to* ~ to in the regime we have studied. 

Here we concentrate on estimating the finite momen- 
tum relaxation rate, 1/rp, of a polaron (gas), due the 
dipole interaction between the two layers, as a function 
of its momentum and of the temperature of the system. 
Such a quantity plays an important role in determin- 
ing the transport properties - in particular the cross re- 
sistance - in double-layer semiconductor structures (see 
e.g., [3l| and reference therein) or double- layer graphene 
systems (e.g., [32J). Indeed in coupled two-dimensional 
electron gases the electron-electron scattering between 
the two layers gives rise to the Coulomb drag effect, where 
a "drag" current is induced purely from the momentum 
exchanges through inter-layer electron-electron scatter- 
ing events. In a different context the polaron momentum 
relaxation rate (or spin-drag [33j]) has already been cal- 
culated for a three-dimensional highly polarized spin-1/2 
Fermi gases, where the interaction is short range 134] . 

In a real polaron experiment one will have in the upper 
layer a gas at very low density ua , much smaller than the 
lower layer density n B . We take the gas in A to have a 
momentum per unit volume = nAmv, with v the 
average velocity of the gas. We also assume that in both 
layers the gases are described by free quasi-particles in 
thermal equilibrium at temperature T whose distribution 
functions are n pB = /[^(e p s-/xs)] and n pA = /[/3(e p A- 
p v - n A )}, with p = l/k B T and f{x) = l/(e x + 1) 
the Fermi distribution. The single particle energies are 
e pA(B) — p 2 /2m. The term p • v boosts the ^4-layer 
distribution function by a velocity v. 

The momentum relaxation rate rp is defined by 



Pa 

dt 



Pa 

TP 



(12) 



We will explicitly give the expressions of Tp for relevant 
cases, closely following the calculation for the Coulomb 
drag between parallel two-dimensional electron gas by 
Jauho and Smith (35| . 

The rate of change of the momentum of the dipolar gas 
in layer A due to the scattering from atoms/molecules in 
the layer B reads 



dP A 

dt 



p,p',q 



«-p-qA)(l - n-p'+qs) - (P P 



q, p' o p' 



q)] ^(epB+ep/A-ep-qA-ep'+qs), 

(13) 
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with w(q) the probability transition for the scattering 
process pa + p' B -> (P q)a + (p' + q)s- Taking the 
continuum limit and introducing the density-density re- 
sponse function for the layer X = A, B 



Eq. (US]) one gets 



d 2 p 



(2tt) 2 e p - e p+q + hu + i0+ 



(14) 



rff 



d 2 q 



qtu(q) 



(2tt)2 

ImxA(q, q • v - w)ImxB(q, w) 



(1 - c^(" _c i' v )(l — e~P") 



,(15) 



Since we are interested in the low temperature behaviour 
we can use the zero temperature response function in 
the previous expression. Moreover for the aim of the 
present paper we just use the golden rule for the collision 
probability with the bare dipolar potential, i.e., w(q) = 
2ir\V(q)\ 2 /h (but see text below). 

Analytical expressions for the relaxation rate are pos- 
sible in some limiting cases. We are interested in partic- 
ular in the dependence of Tp on the various parameters 
at low-T and low-u. 

a. Velocity dependence at T = 0. Assuming zero 
temperature and low velocity such that mv <C kp,A> with 
kp t A the Fermi momentum in the layer A, we can use 
the expression ImxA(B) — m2{jJ / {itPqkp^B^v^k in Eq. 
(1151). We can then write the momentum relaxation time 



as 



1 

TP 



24tt 2 h 



{Xkp^Y \ kp A 



:f,b 



(16) 



b. Temperature dependence at very low-v. In this 
case we find the leading term at small temperature. 
Such a result is usually the one of main interest in solid 
state physics. Assuming that the temperature is low 
enough that the gases in both layers are degenerate, but 
vkF t A < k B T, we can again use the low-frequency ex- 
pression for ImxA(s) and we obtain 



1_ = it (k B T) 2 /V 
T P 32 tfep.B V A 



1 



(AfcF,s) 2 \kF, 



«F,B 



(17) 



General arguments on the phase space restriction for 
the collisions suggest that the leading terms are indeed 
v 2 and T 2 . Within our approximations we find a 1 /A 4 de- 
pendence at finite temperature and at zero-temperature 
a 1/A 6 dependence. 

It is worth noticing that the same finite temperature 
dependence on the layers' distance has been found for 
large enough inter-layer separation for a double layer elec- 
tron gas [35| and graphene [32[ due to the screening of the 
Coulomb potential. Indeed the screened Coulomb poten- 
tial within random phase approximation has the same 
functional form of the dipole inter-layer potential [36| . 



In the random phase approximation the effective dipo- 
lar potential due to (static) screening in the layer- B for 
exchange momentum q < 2kp reads 



VRPA(q) = - 



2nd 2 qe 



-Xq 



1 — r*q 



(18) 



where we used the static value of the polarizability of a 
two dimensional ideal Fermi gas Thus in our case (due 
to the different scaling of the dipolar potential with re- 
spect to the Coulomb one) for large inter-layer separa- 
tion, r*/A < 1, we can consider the unscreened potential 
and find Eqs. ([16]) and ([17]). 



IV. CONSIDERATIONS ON EXPERIMENTS 
WITH ULTRA-COLD TRAPPED GASES 

In the previous sections we consider an uniform two- 
dimensional gas. In most of the experiment on cold gases 
an external parabolic confinement is present. In partic- 
ular a two dimensional configuration is realized provided 
the confinement is strong enough to freeze the motion in 
one direction, e.g., along z. For a cylindrical symmetric 
trapping the potential in the other remaining directions 
x, y can be written as 



V ho (x,y) = ^muj ± (x 2 +y 2 ). 



(19) 



Within Hartree-Fock approximation for the energy of the 
system the density profile of the gas in such a potential 
can be determined by a variational Thomas- Fermi ansatz 
which in 2D is simply an inverted parabola 



n(x,y) 



2N 



x 2 +y 2 



(20) 



where TV is the number of atoms / molecules and R± is the 
variational radius. We notice that since we consider that 
the dipoles are oriented perpendicularly to the plane we 
do not need to consider deformations of the Fermi surface 
(see, e.g., for the 2D case [371. l38j). 

A double layer configuration can also be realized, e.g., 
by means of selecting two wells in an one-dimensional 
optical lattice. For the sake of simplicity in the dscussion 
we assume that the two dimensional confining potential 
is the same for both layers. 

The polaron energy depends on the density of the layer 
B through kp,B- Thus in a local density approach the 
polaron feels an effective external potential 

V p {x,y) = V ho (x,y)+e p (n B (x,y)) (21) 

where we assume the density in the layer B, n B has the 
shape P0|) . 

If it is possible to excite the dipole mode for the polaron 
[39[ | . i.e., the frequency of small oscillations of a dilute 
polaron gas, one will find a value of its frequency lo p 
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larger than the bare value In particular using the 
simple expression Eq. (fTTj) one finds 

I . K\ (, ijgl i M 1/a 

P ~ ^ [ V 2 ei ,, l + | eb |/(2 eF , )J ' 

(22) 

where 6f,o is the value of the Fermi energy at the cen- 
ter of layer B. For small inter-layer interaction we have 
simply lu p /lu± ~ 1 + e^/e^o while it saturates for strong 
interaction. 

In order to measure the frequency Eq. (|22p the dipole 
mode should decay slowly such that enough oscillations 
occur. The damping time of such oscillation is related 
to the momentum relaxation time studied in Sect. V (as 
explained in [IH). In particular if uj p t p 1 the dipole 
mode is collisionless and thus weakly damped, while if 
LOpTp <C 1 the mode is hydrodynamic and therefore over- 
damped. 

It is worth mentioning that for a short-range interact- 
ing spin- 1/2 Fermi gas the frequency of the quadrupole 
(4fj| and dipole [41[ polaron modes and their decay rates 
has been recently measured providing informations on 
the Landau parameters of the unbalanced normal phase 
of the gas 39] and a new insight in the transport prop- 



erties of strongly interacting Fermi system [4l| 



V. CONCLUSION 

In conlusion, we calculated the energy of a dipolar po- 
laron in a double layer configuration for relatively weak 
interactions assuming the bath being in its Fermi liquid 
phase. In the same spirit we determined the momentum 
relaxation rate of a (gas of) polaron as a function of the 
gas' speed and of the temperature of the system. The lat- 
ter one is related to the well known drag effect in a double 
layer electronic structure. Furthermore, we estimated the 
frequency shift and decaying time of the out-of-phase or 
dipolar mode in a trapped configuration, which should 
correspond to a realistic experimental implementation of 
our system. 



Acknowledgments 

Useful discussions with Chiara Menotti, Alexander 
Pikovski and Luis Santos are acknowledged. This work 
has been supported by ERC through the QGBE grant. 



[1] M. A. Baranov, Physics Reports 464, 71 (2008). 

[2] T. Lahaye et al, Rep. Prog. Phys. 72, 126401 (2009). 

[3] K.-K. Ni et al, Science 322, 231 (2008). 

[4] J. Deiglmayr et al, Phys. Rev. Lett. 101, 133004 (2008). 

[5] S. Ospelkaus et al, Science 327, 853 (2010). 

[6] K.-K. Ni et al, Nature 464, 1324 (2010). 

[7] M. Klawunn, A. Pikovski and L. Santos, Phys. Rev. A 
82, 044701 (2010). 

[8] J. R. Armstrong et al., Europhys. Lett. 91, 16001 (2010). 

[9] M. A. Baranov, A. Micheli, S. Ronen, and P. Zoller, Phys. 
Rev. A 83, 043602 (2011). 
[10] A. Pikovski, M. Klawunn, G.V. Shlyapnikov, and L. San- 
tos, Phys. Rev. Lett. 105, 215302 (2010). 
[11] A. C. Potter, E. Berg, D.-W.Wang, B. I. Halperin, and 

E. Demler, Phys. Rev. Lett. 105, 220406 (2010). 
[12] A. Pikovski, M. Klawunn, A. Recati, and L. Santos, 

eprint: arXiv: 1108.5642vl. 
[13] A. Schirotzek, C.-H. Wu, A. Sommer and M. W. Zwier- 
lein, Phys. Rev. Lett. 102, 230402 (2009); S. Nascim- 
bene, N. Navon, K. J. Jiang, L. Tarruell, M. Teichmann, 
J. McKeever, F. Chevy, and C. Salomon, Phys. Rev. Lett. 
103, 170402 (2009); N. Navon, S. Nascimbene, F. Chevy, 
C. Salomon, Science 328, 5979 (2010); B. Frohlich, M. 
Feld, E. Vogt, M. Koschorreck, W. Zwerger, and M. Kohl, 
Phys. Rev. Lett. 106, 105301 (2011). 
[14] F. Chevy and C. Mora, Rep. on Prog, in Phys. 73, 112401 
(2010). 

[15] J. B. McGuire, J. Math. Phys. (N.Y.) 7, 123 (1966). 
[16] S. Giraud and R. Combescot, Phys. Rev. A 79, 043615 
(2009). 

[17] S. Zollner, G. Bruun and C. J. Pethick, Phys. Rev. A 83, 
021603 (2011). 



[18] M. Parish, Phys. Rev. A 83, 051603 (2011). 
[19] M. Klawunn and A. Recati, Phys. Rev. A 84, 033607 
(2011). 

[20] R. Schmidt, T. Enss, V. Pietila and E. Demler, eprint: 

arXiv: 1110.1649. 
[21] R. Combescot, A. Recati, C. Lobo and F. Chevy, Phys. 

Rev. Lett. 98, 180402 (2007). 
[22] S. Pilati and S. Giorgini, Phys. Rev. Lett. 100, 030401 

(2008). 

[23] R. Schmidt and T. Enss, Phys. Rev. A 83, 063620 (2011). 

[24] E. Lipparini, Modern many particle physics (2nd edition), 
World Scientific 2008. 

[25] C.-K- Chan, C. Wu, W.-C. Lee, and S. Das Sarma, Phys. 
Rev. A 81, 023602 (2010). 

[26] For our considerations the influence of the finite width of 
the layers l z <C A is negligible. However, in similiar situ- 
ations it can lead to very interesting physics as discussed 
in the two recent preprints: M. Babadi, and E. Demler, 
eprint: arXiv: 1109.3755v2 and V. Pietila, D. Pekker, Y. 
Nishida, and E. Demler, eprint: arXiv: 1110.0494vl. 

[27] H. A. Bethe and J. Goldstone, Proc. Roy. Soc. A 238, 
551 (1957). 

[28] S. Nagano, K. S. Singwi and S. Ohnishi, Phys. Rev. B 

29, 1209 (1984). 
[29] G. E. Astrakharchik, J. Boronat, I. L. Kurbakov and Yu. 

E. Lozovik, Phys. Rev. Lett. 98, 060405 (2007). 
[30] Stefano Giorgini, private communications. 
[31] A. G. Rojo, J. Phys.:Condens. Matter 11, R31 (1999). 
[32] W.-K. Tse, B. Y.-K. Hu and S. Das Sarma, Phys. Rev. 

B 76, 081401 (R) (2007) 
[33] M. Polini and G. Vignale, Phys. Rev. Lett. 98, 266403 

(2007); R. A. Duine, M. Polini, H. T. C. Stoof, and G. 



7 



Vignale, Phys. Rev. Lett. 104, 220403 (2010). 
[34] G. M. Bruun, A. Recati, C. J. Pethick, H. Smith and 

S. Stringari, Phys. Rev. Lett. 100, 240406 (2008). 
[35] A.-P. Jauho and H. Smith, Phys. Rev. B 47, 4420 (1993). 
[36] T. Ando, A. B. Fowler and F. Stern, Rev. Mod. Phys. 

54, 437 (1982). 

[37] G. M. Bruun and E. Taylor, Phys. Rev. Lett. 101, 245301 
(2008). 

[38] Y. Yamaguchi, T. Sogo, T. Ito and T. Miyakawa, Phys. 
Rev. A 82, 013643 (2010). 



[39] A. Recati, C. Lobo, and S. Stringari, Phys. Rev. A 78, 
023633 (2008). 

[40] S. Nascimbene, N. Navon, K. J. Jiang, L. Tarruell, M. 

Teichmann, J. McKeever, F. Chevy, and C. Salomon, 

Phys. Rev. Lett. 103, 170402 (2009). 
[41] A. Sommer, M. Ku, G. Roati, M. W. Zwierlein, Nature 

472, 201 (2011); A. Sommer, M. Ku, M. W. Zwierlein, 

New J. Phys. 13, 055009 (2011). 



